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Analytical expressions are derived for the melting line and liquid equation of state of normal deuterium near 
the triple point. Melting pressures were measured between the triple point and 20.4 K. These results combined 
with existing pressure measurements along the saturated liquid-vapor curve fix an accurate value, T tp =\ 1 8.723 K, 
for the triple-point temperature. Data for the isothermal compressibility and thermal expansion coefficients of 
the liquid were taken over the temperature and pressure ranges 18.8 to 21.0 K and 4 to 70 bar, respectively. The 
liquid molar volume was measured at nine points below 20-4 K. All liquid PVT data are shown to be internally 
consistent. Measurements of the volume changes on melting are also presented. The heat of fusion and the solid 
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existing data. A critical appraisal is given of all measured thermodynamic quantities in this regime. 
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1. Introduction 

Interest in the properties of the condensed-phase hy- 
drogen systems has continued for well over 50 years. As 
early as 1935, Clusius and Bartholome [l] 1 published the 
first comprehensive thermodynamic study of normal 
deuterium (n-D 2 ) in the triple-point region. These results 
were followed by others [2,3] and included in the 1948 
review article by Woolley et al.[4]. Since that time, there 
have been some additional measurements in this regime 
[5,6], but most efforts have been toward establishing the 
pressure-volume-temperature (PVT) relationship over 
extended ranges of temperatures and pressures. 

By 1959, an incomplete collection of fluid-phase den- 
sity data was available at various intervals in the 

About the Authors, Paper: The work reported on 
was performed at the Los Alamos National Labora- 
tory where L. A. Schwalbe, a physicist, remains and 
from which E. R. Grilly, also a physicist, is retired. 
The work was supported by the Fusion Target Fabri- 
cation Group at Los Alamos. 

'Figures in brackets indicate literature references at the end of this 
paper. 



temperature range from the triple point to 420 K and at 
pressures to 3 kbar. Prydz [7,8] critically reviewed these 
results, fit them to a modified Strobridge equation, and 
derived analytical expressions for selected thermo- 
dynamic properties. 

More recently Mills et al. [9] measured simulta- 
neously the molar volume and ultrasound velocity of 
fluid n-D 2 in a piston-cylinder apparatus. The data 
spanned temperatures of 75 to 300 K and pressures of 2 
to 20 kbar. The results were fit to a Benedict-type equa- 
tion of state, and analytic forms were derived for both 
constant-volume and constant-pressure heat capacities 
over these ranges. In a companion to this work, Lieb- 
enberg et al. [10] reported measurements of the melting 
line from 4 to 19 kbar and the corresponding changes in 
the molar volume and longitudinal sound velocity. 
From these data, they derived the entropy change on 
melting, the adiabatic compressibility for the solid and 
other thermodynamic properties. 

Most work on solid deuterium has been done at low 
temperatures and low pressures although the volume- 
pressure dependence at 4 K was measured to 20 [1 1] and 
25 kbar [12]. Recently, solid-phase isochoric equation- 
of-state data were taken to 2 kbar [13,14]. 
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The purpose of this study is to supplement existing 
PVT and melting data with accurate and precise 
measurements 2 in the region of the phase diagram near 
the triple point. Data are presented for the pressures and 
liquid molar volume of n-D 2 along the melting line from 
the triple point to 20.4 K. The isothermal com- 
pressibility and thermal expansion coefficients of the 
liquid are given over the same temperature range and at 
pressures between 4 and 70 bar. The combined liquid 
data are fitted to an empirical equation of state and are 
shown to be internally consistent. 

Measurements of the volume changes on melting are 
also presented. The heat of fusion and the solid molar 
volume at melting are deduced from these data. 
Throughout the discussion, we have included detailed 
comparisons of our results with existing data and have 
attempted to provide a critical appraisal of all measured 
quantities in this regime. 

2. Experimental 

Our studies on n-D 2 were carried out with essentially 
the same apparatus as that used earlier for similar work 
[15-17] on 3 He and 4 He. The pressure cell, which was 
referred to in previous descriptions [16,17] as "cell I," 
consisted of three BeCu diaphragms that were welded 
together at their circumference and separated by 0.3 mm 
gaps (see fig. 1). The lower gap was connected to a 
capillary tube and remained open to a room-temperature 
helium gas-handling system. With this arrangement P u 
the pressure of the lower-cell volume, could be adjusted 
or measured directly. The upper cell gap served as the 
deuterium sample chamber; a fixed quantity of material 

2 Full data are available from the authors. 



could be maintained in its volume V u by sealing the 
low-temperature valve. The sample pressure P a was de- 
termined from the deflection of the upper diaphragm. 
Originally the deflection measurement was made with a 
linear differential transformer [15], but the capacitance 
technique was adopted because of its greater sensitivity 
[17]. 

Changes in the volume of the sample chamber depend 
upon changes in P m P,, and the bath pressure P b accord- 
ing to the equation 



A F»=(S 11 +5 l )AP u -SM-5 B AJ» b 



(1) 



where iS u and S\ are pressure sensitivity factors of the 
volume changes to changes in the upper and middle 
diaphragm displacements, respectively. 

The pressure sensitivity factors for the cell were not 
constants but varied reproducibly by 1 or 2% as P { and 
P a varied between and 70 bars, the highest working 
pressures in this experiment. For the earlier helium mea- 
surements, these variations did not critically affect the 
results, but for the present work it was first necessary to 
calibrate the system more precisely. To do this, three 
separate sets of calibration measurements were run. In 
each, the procedure first involved loading the sample 
chamber with liquid p-H 2 at T= 20.00 K and then 
"bleeding" small amounts from the cell through the 
valve and into a calibrated standard volume held nomi- 
nally at room temperature. The quantity of material 
drawn from the cell was then measured and the corre- 
sponding changes in the sample volume AV V were calcu- 
lated using the p-H 2 density values tabulated by Good- 
win et al. [18]. According to eq (1), isothermal changes 
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in V u under the constraints AP u =0, AP a =AP u and 
APj=0 yield the factors S t , S a , and S u +S\, respectively. 

A consistent set of pressure sensitivity parameters was 
obtained at r=20.000 K for the entire working range of 
both upper- and lower-cell pressures. The quantities S a 
and S\ that were used in the data analysis were accurate 
to ±0.1%. Extrapolations of total volume mea- 
surements to P u =P]=P b =0 absolute pressure deter- 
mined Puo, the zero-pressure volume of the sample 
chamber. The latter calibration was made at two sepa- 
rate times; initially, before any of the deuterium work 
was begun and again near its completion about 18 
months later. The first set of measurements obtained an 
average value 0.37477 cm 3 for six separate runs; for the 
second set the average of five runs was 0.37426 cm 3 . In 
the data analysis we used the mean of these two values, 
^0=0.37452 cm 3 , and set an accuracy limit of ±0.15% 
for this quantity. This error includes the ±0.1% uncer- 
tainty of the original p-H 2 density data [18]. 

Equation (1) has no explicit temperature dependence 
included. In this range, there are negligible effects intro- 
duced by the temperature coefficient of the cell mate- 
rial. The thermal expansion of BeCu should be compara- 
ble to that of Cu. Rubin et al. [19] reported a linear 
expansion coefficient of 6X 10 -6 K _1 for copper at 25 K. 
The effects introduced by thermal expansion are, there- 
fore, about an order of magnitude smaller than the 
quoted uncertainties in our data. To check this explic- 
itly, we measured thermal expansion coefficients of liq- 
uid hydrogen in our cell. Hydrogen data were taken at 
0.4 K intervals between 19.0 K and 21.0 K along six 
isobars (5.48, 14.76, 28.52, 42.28, 56.05, and 69.81 bar). 
The measured values agree with those that Roder et al. 
[20] derived from the original density data of Goodwin 
et al. [18]. 

The uniform sample temperatures required for this 
study were maintained by holding the pressure cell as- 
sembly immersed in a bath of liquid 20 K-equilibrium 
hydrogen (e-Hi). Temperatures were controlled by reg- 
ulating the bath pressure with a Baratron 3 Type 
77H-1000 metal diaphragm gauge (0-1000 Torr range). 
Bath pressures were obtained directly with an 
18-mm-bore mercury manometer. The mercury column 
heights were measured with a cathetometer, and stan- 
dard corrections for the temperature and local gravity 
constant were applied to these readings. From these 
data the sample temperatures were calculated with the 



saturated liquid-vapor pressure equation of Souers et al. 
[21]: 
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In ft(torr)= 10.5741 1-101.3378 T~ l 
+ 5.432005 XlO- 2 T-\. 105632X10"- 4 T 2 . 



(2) 



This function is identical to that defined by the Inter- 
national Practical Temperature Scale of 1968 (IPTS-68) 
[22] although it is expressed in somewhat different form. 
Temperature precisions of slightly better than 1 mK 
were attainable with this procedure. From run to run, 
with the possible exceptions noted below, the scale 
seems to have been reproducible to within about ±2 
mK. 

3. Results and Discussion 
3.1 Melting Pressures 

The melting pressures of n-D 2 were measured by clos- 
ing the cell valve on a liquid-phase sample and in- 
creasing Pi until the capacitance readings showed the 
sample to be in the two-phase region. At each tem- 
perature, capacitance readings were recorded for 
several different P t values. Finally, the valve was opened 
and the capacitance-pressure calibration was rechecked 
against the piston gauge. The results are given in table 1. 
Measurements were taken during four runs, each lasting 
one week; each run represents a separate filling of e-Hj 
refrigerant 

It is important to note that all melting pressures re- 
ported here correspond to the "first-freeze" character- 
istics of the sample. The pressure-temperature mea- 
surements were recorded for small solid fractions 
(<5%) of the two-phase mixture. In a later section, we 
discuss in more detail the problems involved with pre- 
cision measurements of the sample pressure when larger 
relative quantities of solid are present. 

Values for the slope of the melting curve are plotted 
in figure 2. Those obtained from run numbers 1, 2, 3b, 
and 4 were derived from consecutive melting pressure 
measurements. Those labeled "3a" were measured di- 
rectly by temperature cycling a single sealed sample. 
The plot indicates a nearly linear temperature de- 
pendence for the slope. The melting line P m (T m ) should, 
therefore, be adequately represented by a quadratic 
function over this limited range. For the least-squares 
fit, we chose an equation of the form 



pjr m )=p lp +A l {T m -T ip ) +A 2 (T m - jy 1 



(3) 



where A\ and A 2 are fitting constants and P^ and T tv are 
parameters that represent the pressure and temperature 
values at the triple point. Our best fit gave >li= 38.884 
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Table 1. Measured pressure and temperatures of normal deuterium along the melting line. 



Hun No. 1 r m (K) P m (bar) 



Run No. 2 T m (K) /> m (bar) 



Run No. 3 r m (K) P m (bar) 



Run No. 4 T m (K) P m (bar) 



20.407, 


69.61 


20.206o 


60.75 


20.006s 


52.40 


19.805, 


44.02 


19.605, 


36.10 


19.404a 


27.88 


19.198 2 


19.61 


20.306, 


65.29 


20.107, 


56.83 


19.906o 


48.12 


19.704 6 


39.79 


20.408, 


69.87 



19.504i 


31.90 


19.302* 


23.73 


19.1022 


15.72 


18.902o 


7.81 


20.007 


52.57 


19.804 2 


44.14 


19.603 3 


35.93 


19.402, 


27.75 


20.407 4 


69.44 


20.305( 


65.10 


20.205 s 


60.83 


20.105+ 


56.63 


19.906 2 


48.05 



19.703o 


40.03 


20.407 


69.44 


20.204j 


60.86 


20.004, 


52.48 


19.802 2 


44.08 


19.192, 


19.40 


19.000, 


11.75 



20.306, 


65.16 


20.103, 


56.62 


19.902, 


48.24 


19.699 s 


39.86 


19.504, 


31.88 


19.301, 


23.68 


19.100, 


15.64 


18.8992 


7.69 


19.404, 


27.79 


19.195! 


19.38 


19.000, 


11.65 


18.829i 


4.88 


19.603 8 


36.04 


18.828, 


4.96 



bar/K, ^,= 1.078 bar/K 2 , and T^= 18.7067 K with 



/> tp =0.17 bar. 



Grilly [23] reported a triple-point temperature of 
18.73 K. The result obtained here 18.7067 K is consid- 
erably lower. We attempted to fit eq (3) by varying only 
Ai and^ 2 with T tp = 18.73 K and ^=0.17 bar held fixed, 
but the results showed large systematic deviations. We 
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Figure 2— The slope of the melting curve. Data from runs 1, 2, 3b, and 
4 were obtained by subtracting consecutive measurements of the 
melting pressure. Data from run 3a were measured by temperature 
cycling a single sealed sample. The solid line was obtained by 
evaluating eq (4). The dashed line is the slope of the melting func- 
tion given by Clusius and Bartholome [1]. 



believe the difference is due to impurity effects in the 
actual fixed-point values of our sample lot. Mass spec- 
trometric analysis of the sample material showed 0.75% 
HD present with no detectable traces of H 2) N 2 , or H 2 
at the 50 ppm sensitivity level. Sample impurities in low 
concentrations affect both triple-point parameters, but 
within the resolution limits of this experiment, the tem- 
perature T tp is much more sensitively dependent on im- 
purity effects than is the pressure Ftp- Therefore, in the 
final calculation we retained the fixed value i> tp =0.17 
bar and allowed only T lv to vary along with A\ and A 2 . 

Equation (3) with the above parameters describes the 
melting line for our sample, but for pure n-D 2 , a pressure 
correction must be applied. If this correction is propor- 
tional to the HD concentration and to the difference in 
pressure between the HD and D 2 melting curves, 0.6 bar 
should be subtracted from P m in eq (3) to obtain the 
melting pressure of pure n-D 2 . 

A plot of the deviations of our data from eq (3) is 
shown in figure 3. With the exception of a few points, 
mainly all from the first run, the data appear to lie within 
±0.1 bar of the smooth curve. The precision of our 
measurements is comparable to that of Bereznyak and 
Sheinina [24], who measured the melting lines of e-H 2 , 
n-H 2 , and e-D 2 over roughly the same pressure interval. 
Their P m measurements on e-D 2 exceed ours on n-D 2 by 
1.5 bar at 18.8 K and 2.3 bar at 21.0 K. (This difference 
is similar to that observed between the melting lines of 
e-H 2 and n-H 2 .) However, they observed a slope discon- 
tinuity in their e-D 2 melting curve at T= 19.02 K and 
i> = 14.2 bar. Our data show no anomalous behavior 
anywhere in the range 18.7067 to 20.400 K. 

An analytic form for the slope of the melting curve is 
obtained by differentiating eq (3) with respect to T m : 



^=^,+2^ 2 (r m -r tp ). 



(4). 



230 



0.4 



0.2 



E 
Q. 



0.0 



-0.2 - 



-0.4 - 



1 1 1 1 r 

NORMAL DEUTERIUM 
MELTING PRESSURE 
DEVIATIONS 



'9 * 



A V 

o 



o O 



O RUN No. I 
O RUN No. 2 
A RUN No. 3 
V RUN No. 4 

__| [__ 



19 



20 



TEMPERATURE (K) 



Figure 3 — The deviations of the melting pressure data from the 
smooth curve defined by eq (3). 

We have included in figure 2 a plot of eq (4) with the 
values A U A 2 , and 7^ found above. This function repre- 
sents the data for the slope of the melting curve to 
within ±0.5%. 

Other melting data are available for comparison. 
Among the earliest of these are the Clusius and Bar- 
tholome [1] results on n-D 2 . These measurements extend 
from the triple point up in pressure to about 100 bar. 
Clusius and Bartholome employed a technique in which 
the liquid sample was held at nominally constant pres- 
sures while the surrounding n-H 2 bath pressure was 
gradually lowered. The freezing transition temperature 
was measured by observing the effect of the latent heat 
of the sample on the cryogen bath pressure. Their 
quoted temperatures were obtained from the measured 
cryogen vapor pressure and the centigrade-scale, 6(P), 
vapor-pressure equation of Keesom et a]' [25]. Absolute 
temperatures, according to Bartholome [2], were then 
calculated from the expression T(K)= 0+273. 15. Clus- 
ius and Bartholome fit their results to a Simon-type 
melting equation. This curve, corrected to the IPTS-68, 
is seen to lie about 2 bar higher than the present mea- 
surements. 

Later studies have considerably extended the pressure 
range for the deuterium melting line. Chester and Dug- 
dale [26] measured n-H 2 and n-D 2 to 2.8 kbar, but their 
results cannot be used for direct comparison here be- 
cause the data were expressed only as melting pressure 
differences between the two isotopes at the same tem- 
perature. They observed a nearly constant difference of 
167 ±8 bar between 25 and 57 K. 



Further results were reported by Mills and Grilly 
[27,28] who used a blocked-capillary technique. Their 
first data set [27], on n-H 2 and n-D 2 only, was limited in 
pressure to 1.5 kbar because of embrittlement and crack- 
ing failures in the measurement apparatus. However, a 
later, improved system [28] allowed them to measure 
pressures of n-H 2 , n-D 2 , and n-T 2 to 3.5 kbar. Their 
deuterium samples had about the same HD concen- 
tration as ours. Below 70 bar, the limit of the present 
study, they obtained three points which we corrected to 
IPTS-68 by adding 6 mK to their temperatures. Com- 
parison of these results with eq (3) shows their P m mea- 
surements falling lower than our smooth curve by 0.4 to 
0.8 bar. A comparison of the melting curves obtained for 
n-H 2 and n-D 2 shows a difference of about 170 bar near 
20 K. 

Goodwin [29] proposed the empirical function 

(P m -P, p y(T m -T iv )=A exp(-a/T m )+BT m +C (5) 

as a universal equation for the melting lines of the three 
hydrogen isotopes. For C =0, the constants A, a, and B 
were fixed from existing normal and para -hydrogen 
melting data. The equation was then applied to n-D 2 and 
n-T 2 using their respective triple-point parameters, 7^ 
and P lp . Goodwin adjusted the values of C to obtain best 
fits to the melting data of the two heavy isotopes. At 
pressures exceeding 100 bar, the function for deuterium 
was made to agree with the melting equation published 
by Mills and Grilly [28] to better than 0.5%. Goodwin's 
equations for n-H 2 and n-D 2 differ in pressure by 169 bar 
in the 20 K region. 

Several years after this work was done, Younglove 
[30] obtained additional melting data on p-H 2 and found 
these to be systematically offset from Goodwin's origi- 
nal function [29,31]. At low pressures, Younglove's 
melting pressures are about 0.4 bar lower than Good- 
win's, and at high pressures (ca. 130 bar) they are 1.2 bar 
lower. It is, therefore, interesting to compare our mea- 
surements on n-D 2 with eq (5). In doing so, we find the 
temperature-corrected [32] form of Goodwin's equation 
to be about 0.5 to 0.8 bar lower than eq (3). However, 
the earlier function refers to a triple-point temperature 
of 18.72 K, which more nearly represents the pure n-D 2 
value than does our direct measurement of 18.7067 K. 
The impurity correction of —0.6 bar discussed above 
brings the two melting curves into nearly complete cor- 
respondence. We suggest that Goodwin's equation de- 
scribes the melting pressures of pure n-D 2 to about ±0.2 
bar at pressures below 70 bar. 

The only other melting curve available is that of Lieb- 
enberg et al. [10]. Although the range of validity of this 
experimental curve extends to 20 kbar, there are large 
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(=40 bar) discrepancies in the low-pressure range of 
present interest. 

3,2 Triple-point Temperature 

As we have mentioned, the physical significance of 
r :p in eq (3) is that it represents an empirical triple-point 
temperature For our sample material. A review and dis- 
cussion of existing triple-point data for deuterium is 
presented separately [33]. 

Most of the earlier triple-point measurements were 
derived from saturated solid- and liquid-vapor pressure 
data. With this approach, separate empirical functions 
are fitted to the data in both solid and liquid temperature 
regimes. Simultaneous solution of these functions then 
defines the two triple-point parameters. In general, this 
technique allows P 1p to be fixed with precisions com- 
parable to those of the direct pressure measurements. 
However, the temperature value at the intersection 
point is difficult to find closer than about 20 or 30 mK 
[33]. The greater relative uncertainty in 7^ results from 
the small difference in slope between the solid and liquid 
vapor pressure functions. 

Greater precision for T lp can be obtained with supple- 
mentary melting data because the slope of the melting 
curve is very much larger than that of either the solid or 
liquid vapor pressure functions. Even with a relatively 
approximate estimate of the triple-point pressure, the 
temperature T, p should be determined quite precisely. 
(For example, an error of 25% in P tp would offset 7^ by 
only about 0.001 K.) By similar argument it is clear that 
imprecise estimates of T lp taken as fixed parameters in eq 
(3) should introduce systematic deviations in the re- 
sidual plot corresponding to figure 3. We examined this 
effect explicitly and found that significant deviations are 
apparent in the fitting results for all fixed T ip values that 
lie outside the range 18.7067± 0.002 K. 

This value cannot be taken as the triple-point tem- 
perature of pure n-D 2 because the HD impurity intro- 
duces some offset. Bereznyak et al. [34] published melt- 
ing diagrams for mixtures of p-H 2 /o-D 2 , p-H 2 /HD, and 
o-D 2 /HD. The melting -temperature function of isotopic 
concentration for the o-Dj/HD system was found to be 
nearly linear between the triple points of the two pure 
elements. It is, therefore, reasonable to assume that a 
small concentration of HD-impurity e(HD) in other- 
wise pure n-D ; will shift the observed triple^point tem- 
perature from the pure-component value by an amount 



Ar tp =c(HD)[r tp (HD)- r lp (n-D 2 )]. 



(6) 



Substituting 7, P (HD)= 16.60 K, T ta (n-D 2 )= 18.73 K, and 



c(HD)=0.0075 into eq (6) yields a temperature shift of 
A7* tp = -0.016 K. 

There is another possible source of systematic error. 
We noted earlier that our melting data are actually ob- 
servations of the "first-freeze" behavior of the material. 
Grenier and White [35] similarly measured "first- 
freeze" temperatures of deuterium. But, they also mea- 
sured melting temperatures of samples that contained 
larger relative fractions of the solid phase, and for these 
they reported differences of "a few millikelvin." The 
effect was attributed to variations of the ortho-para com- 
positions of the solid and liquid phases. While this sug- 
gestion may be true, we feel the data are insufficient to 
warrant any firm conclusions. Recall that our direct 
measurements of the melting curve slope (labeled "3a" 
in Fig. 2) gave results that are consistent with those 
derived from the differences of consecutive "first- 
freeze" melting pressures. In our experiment, the melt- 
ing characteristics were affected little by the relative 
phase composition of the sample. 

In view of these uncertainties, we feel that a triple- 
point temperature of 18.723 ±0.010 IC for pure n-D 2 is 
the most reasonable estimate from our measurements. If 
this value is used in eq (3), the calculated P m values will 
have been corrected for the HD impurity. 

3.3 Liquid Isothermal Compressibility 

The isothermal compressibility is calculated from 
the measured change in the sample pressure AP H that 
results when the lower chamber pressure is changed by 
an amount LP\ at constant temperature. The expression 
used in this calculation 

'-KiMM'-g)] (7 > 

is easily derived from eq (1). In this equation, V a repre- 
sents the average sample volume for the measurement 
For small changes in P t and P& the sensitivity factors S a 
and S, may be taken directly from the calibration curves. 
However, to obtain more precise compressibilities for 
the liquid, we chose to measure over larger pressure 
intervals. This procedure requires the sensitivity param- 
eters in eq (7) to consist of rather complicated weighted 
averages of the measured S u and Si functions. 

A plot of the liquid isothermal compressibility mea- 
surements along three isotherms is shown in figure 4. 
The values obtained by extrapolating eight isotherms to 
the melting curve are represented by circles in the plot. 
The estimated accuracy is about ±1%. 

The functional behavior of these data suggests that 
we consider a liquid equation of state of the form: 



232 



10.0 



6.0 



LIQUID NORMAL DEUTERIUM 
ISOTHERMAL COMPRESSIBILITY 



T = 21.00 K 




T = 20.00 K 



* INCREASING PRESSURE 

T DECREASING PRESSURE 

o ISOTHERM EXTRAPOLATED VALUE 
AT THE MELTING PRESSURE 

' I I | l__ 



10 



20 30 40 50 

PRESSURE (bars) 



60 



70 



Figure 4— Isothermal compressibility data for liquid n-Di. Solid 
curves are isotherms calculated from eq (9) and the fitting constants 
listed in table 1. The dashed curve represents eq (9) evaluated along 
the melting line. 



V(T,P)=V (T)[P+P (T)]-° (8a) 

Vo(T)=exp[b (i +b l (,T-T ip nb 1 (T-T tr> y 

+6 3 (r-r tp ) 3 ] (8b) 

^o(r)=co+c 1 (r-r, p )+ C2 (r-7 , lp ) 2 (8c) 

where 7 T tp = 18.7067 K, and a and the subscripted letters 
b and c are all constants to be determined by fitting eq (8) 
and their derivative functions to the experimental data. 
For the isothermal compressibility, we obtain the fol- 
lowing relationship: 



*?*>--{*&),- 



P+PW 



(9) 



A fit of eq (9) to the compressibility data gives the 
numerical values for the constants a, c 0l c u and c 2 listed 
in table 2. In figure 4, the smooth curves depicted for 
three isotherms and the liquid compressibility along the 
melting line were calculated from eq (9) and these con- 
stant values. 

To our knowledge, there are no published com- 
pressibility data in this temperature and pressure range 



Table 2. Coefficients to eqs (8), (9), and (11) that describe the equation 
of state for the liquid. 





a =0.12 


c = 152.0 




ba= 3.7408 


c, = -13.2 




6i = 0.OO17 


c 2 =0.50 




bi = 0.00028 






63=0.00001 







that could be used for direct comparison. But, there 
have been several measurements of liquid molar volume 
at pressure and some estimates of compressibilities can 
be calculated from these by taking first differences. Bar- 
tholome [2] published volume data along three iso- 
therms between 19.5 K and 21.0 K. Compressibilities 
obtained from Bartholomew raw data are too scattered 
to allow any meaningful comparison, but those calcu- 
lated from his smooth curves show qualitatively correct 
pressure and temperature dependences. Quantitatively, 
the agreement with our results is less satisfactory; com- 
pressibilities derived from the earlier work typically lie 
10 to 14% higher than those calculated from eq (9). 

Friedman et al.[6] published molar volume data for 
liquid normal deuterium at seven temperatures between 
20.3 and 38.0 K and at pressures to 100 bar. Their vol- 
ume measurements along the 20.33 K and 20.31 K iso- 
therms were shown to agree reasonably well with the 
earlier but less precise Bartholome [2] data. Com- 
pressibilities derived from the Friedman 20.33 K data 
are in ± 1 % agreement with values found from eq (9). 
Those calculated from their isotherms at 20.31 K. and 
23.52 K. scatter about our fitted smooth curves with an 
average absolute deviation of about 6% and 4%, re- 
spectively. 

3.4 Liquid Thermal Expansivity 

The isobaric thermal expansion coefficient a is calcu- 
lated from the measured change in the lower-cell pres- 
sure Ai>[ that is required to maintain a constant sample 
pressure P u when the temperature of the system is 
changed by an amount AT. From eq (1) it follows that 



a= 



V\dT} P V a AT v\dTj U; 



where V a is again the average sample volume for the 
measurement, and dP b /dT'v& the temperature derivative 
of the liquid p-H 2 vapor pressure calculated from eq (2). 
Figure 5 is a plot of the liquid expansion data obtained 
at six different pressures. We estimate the experimental 
uncertainty for these measurements to be about ± 1 or 
2%. The data precision was limited to some extent by 
the temperature measurements but mostly by a slight 
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Figure S — Thermal expansion data for liquid n-D 2 . Solid curves were 
calculated from eq (11) and the fitting constants listed in table 1. 
The dashed curve represents eq (1 1) evaluated at saturated liquid- 
vapor pressures and melting as indicated. Open circles are calcu- 
lated values from the fitted equation by Kerr [5]. 

nonlinearity in the pressure readout of the TI-144 bour- 
don gauge used for measuring P h Occasionally some 
hysteresis effects were observed in the deflection of the 
cell diaphragm although these only occurred for the 
initial points in either warming or cooling series of mea- 
surements. Errors of about 6% were observed for such 
artifacts, and these few data were discarded. 

Differentiation of eq (8) yields the following expres- 
sion for the thermal expansivity: 

a(7YP)=&, +2^(7- T tv )+3b 3 (T- T tp f 

-/?(7\.P)[c l + 2c 2 (7'-T' tp )] (11) 

where 7^=18.7067 and P(T,P) is calculated from eq 
(9). The constants a and the subscripted c's are those 
determined from the compressibility measurements. 
These constraints assure that the identity 

£),--(£), 

holds and, thereby, provide a condition of internal con- 



sistency that is automatically met if the expansivities can 
be fitted to eq (11) by adjusting only b\, b 2 and b y 

Equation (11) was fitted to the thermal expansivity 
data to fix the values in table 2 for the three additional 
constants. Smooth curves corresponding to the pres- 
sures of the six experimental isobars were calculated 
from eq (11) and are plotted with the data. The dashed 
lines in figure 5 represent the values along the vapor- 
pressure and melting curves as indicated. 

As before with the compressibilities, we found our 
thermal expansivity measurements to be the only ex- 
plicit data in this range. Here again, we have derived 
estimates from the smooth curves that were fitted to 
independent measurements of the molar volume [2,5]. 

Kerr [5] measured the molar volume of the liquid 
along the vapor pressure curve. He combined his results 
with those of Clusius and Bartholome f 1] and fitted them 
to a quadratic function of temperature. We calculated 
thermal expansion coefficients from this function by dif- 
ferentiation. The slight correction to constant pressure 
was made with eq (9) and Grilly's [23] vapor pressure 
function. The results of this calculation, included in fig- 
ure 5 as the series of open circles, agree reasonably well 
with eq (11) near 19 K but deviate systematically with 
temperature to more than 2% at 20.5 K. 

Bartholome [2] presented measurements of the liquid 
molar volume as a function of pressure along three iso- 
therms. Again, we consider the smooth curves that he fit 
to his raw data. The differences between these isotherms 
at constant pressure were taken to calculate average 
thermal expansivities in the two temperature intervals. 
At low pressures there is satisfactory agreement with 
our data. (Here, the consistency is reasonable if it is 
recalled that Bartholome calibrated the volume of his 
pressure cell along the vapor curve against the molar 
volume data of Clusius and Bartholome. The latter mea- 
surements were included by Kerr [5] in his quadratic fit, 
and thermal expansion coefficients derived from this 
function were shown to correspond closely to our 
extrapolated measurements.) However, with increasing 
pressure, the average thermal expansion coefficients cal- 
culated from Bartholome's data decrease much more 
rapidly than ours do. At 70 bars they fall to about 50% 
of our results. 

3.5 Liquid Molar Volume 

Data for V\ m , the liquid molar volume along the melt- 
ing curve, were obtained by a two-step operation. First, 
the volume of a sealed liquid sample was measured at a 
given temperature and two-phase pressure P a . The sam- 
ple quantity was then determined by expanding the gas 
into a standard volume at room temperature. The latter 
measurement followed essentially the same procedure 
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as that used for the calibration measurements of the cell 
volume and pressure-sensitivity factors. 

The results are listed in table 3. We estimate their 
accuracy and precision to be about ±0.035 and ±0.003 
cmVmole, respectively. For all of the measurements 
except one, the data were taken within about 6 h from 
the time of the sample loading. The point at !T m = 19.000 g 
K was obtained from a sample that had been loaded and 
held at that temperature for approximately 17 h before 
the measurement. 

The results are plotted together as a function of tem- 
perature in figure 6. Also included in the figure are three 
points that Woolley et al. [4] had compiled from earlier 
work [1,2]. The liquid molar volume measurements of 

Table 3. Measured liquid molar volumes along the* melting curve. 
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Figure 6— Measurements of the molar volume of liquid n-D 2 at melt- 
ing. The solid curve was calculated from eq (8) at the temperatures ■ 
and pressures prescribed by eq (3). A ±0.15% error is indicated for 
these results. The liquid density function of Prydz [7] was evaluated 
at the same points and the results are represented in the plot by the 
dashed curve. Also included are three points that were tabulated by 
Woolley et al. [4], 



Bartholome were discussed in the two previous sec- 
tions. To obtain the points at 19.70 and 20.3 1 K, Woolley 
evaluated Bartholomew smoothed functions at the tem- 
peratures and pressures prescribed by the melting curve 
of Clusius and Bartholome [1], The triple-point value, 
23.14 cm Vmole, was taken directly from the work of 
Clusius and Bartholome [1], 

Prydz [7] used these and other earlier data to Fit an 
equation of state for liquid normal deuterium. We evalu- 
ated the Prydz density function at temperatures and 
pressures corresponding to our observed melting line. 
The resulting curve, which is also plotted in figure 6, 
adequately represents the earlier data. However, our 
measurements indicate that the liquid compression 
along the melting curve is not nearly as great as the 
Prydz function would suggest. 

The liquid molar volume data were taken with the 
parameters a, b u b 2 , 6 3 , c 0> c u and c 2 in table 2 to fix the 
remaining constant b . With the full set of fitting param- 
eters, eq (8) can be used to calculate liquid n-D 2 molar 
volumes for any temperature and pressure within or 
even somewhat beyond the ranges of the measurements. 
The curve depicted in Figure 6 was calculated along the 
melting line defined by eq (3). The average absolute 
deviation of the measured data from the calculated 
curve is 0.0026 cm'/mole, which is comparable to the 
estimated precision of the measurements. 

We have shown in this and in the two preceding 
sections that the isothermal compressibility, thermal ex- 
pansion, and molar volume data for the liquid are all 
well described by a single empirical equation of state. 
This is a sufficient demonstration of the internal consis- 
tency of these measurements. It remains now to com- 
pare the predictions of this equation with independent 
data, but First some comment should be given about the 
effects of sample contamination. 

Our volume data and their representation, eq (8), 
Were presented without any correction for impurity ef- 
fects. The choice is somewhat arbitrary as to whether or 
not these should be formally included. Some of the exist- 
ing data were given as corrected values [5,36], but oth- 
ers [1,2,6] were published without corrections as we 
have done. The effects are not entirely negligible. For 
our case, the only significant impurity was the 0.75% 
HD discussed above. If we assume the volume cor- 
rection is proportional to the HD concentration and to 
the difference between the molar volumes of HD and 
D 2 , we have 



AK 1= c(HD) [K 1 (HD)~F l (n-D 2 )]. 



(13) 



The most recent density measurements on HD are those 
of Rudenko and Slyusar' [37], which extend from 16.6 to 
35.5 K along the vapor curve. These data, combined 
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with extrapolated values from eq (8), yield a volume 
difference of about 2 cmVmole. The correction to eq (8) 
is then toward smaller values by an amount 0.015 
cmVmole. With the possible exception of some of the 
earliest work [1,2], the extent of impurity effects in pre- 
vious data is comparable to that estimated for our own. 

Bartholome [2] published the earliest measurements 
of the molar volume of n-D 2 at pressure. Bartholomew 
raw data scatter about the fitted smooth curves with an 
average absolute deviation of 0.025 cmVmole. The pre- 
cision of these measurements is, therefore, roughly an 
order of magnitude less than that typical of more recent 
work [5,6]. The average absolute deviation of Bar- 
tholome's data from eq (8) is 0.073, 0.035 and 0.058 
cmVmole for the isotherms at 19.723, 20.346 and 21.032 
K (IPTS-68), respectively. Although these differences 
are probably no larger than the combined error for the 
two experiments, the earlier data show significant sys- 
tematic deviations from eq (8). As has been mentioned 
the compressibilities derived from Bartholomew work 
are larger than our measurements by 10 and 14% at 
:T=20.346 and 21.032 K (IPTS-68), respectively. The 
discrepancies may be the result of sample contamination 
which we shall discuss later in more detail. 

The later measurements by Friedman et al. [6] pro- 
vide us with a more critical basis for comparison. Of 
particular interest here are their tabulated results along 
the isotherms at 20. 3 1 and 20. 3 3 K. We plotted these two 
sets of volume data as functions of pressure and found 
them to agree nearly identically. To compare these with 
our results, we next included the corresponding iso- 
therm calculated from eq (8). The plot of the calculated 
values parallels the Friedman measurements but is uni- 
formly lower in magnitude by about 0.05 cmVmole. 
This discrepancy is larger than our estimated error but 
probably smaller than the combined experimental un- 
certainties. Although Friedman et al. gave no explicit 
error limits, an accuracy comparable to our own 
(±0.035 cmVmole) is reasonable. With this assumption 
the combined error is ±0.07 cmVmole, and their vol- 
ume measurements at 20.3 K are shown to agree with 
our equation of state. A similar comparison was made 
with their data at 23.52 K. At this higher temperature, 
the difference between the earlier measurements and the 
values calculated from eq (8) is no greater than 0.035 
cmVmole throughout the entire pressure range. 

The Friedman data were presented without any ex- 
plicit impurity corrections although the authors gave an 
impurity level of 0.2 at % hydrogen (equivalent to 0.4% 
HD), roughly half the HD-concentration of our sam- 
ples. This result seems improbable because most avail- 
able deuterium has the same impurity as our sample, and 
they apparently made no special effort to purify their 



material. Ultimately, the accuracy of their impurity 
analysis (whether they had 0.4 or 0.8% HD) does not 
change the overall conclusion. The correspondence be- 
tween their data and ours is not significantly affected by 
differences in contamination levels. 

In a previous section, we mentioned that independent 
measurements of the liquid volume are available along 
the vapor pressure curve. These should provide useful 
comparisons, but because our lowest working pressures 
were between 4 and 5 bar, it is necessary to extrapolate 
our results somewhat. To do this, we calculated a molar- 
volume function of temperature by combining Grilly's 
[23] vapor-pressure equation for liquid n-D 2 with eq (8). 
A plot of this curve is shown in figure 7. Also included 
are the values Clusius and Bartholome [1] reported in 
the range 18.80 to 20.53 K, the measurements Kerr [5] 
made between 19.5 and 24.5 K, and the values found by 
extrapolating the data of Bartholome [2] to the vapor 
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Figure 7— Existing molar volume data for liquid n-D 2 along the satur- 
ated vapor pressure curve. Direct measurements by Clusius and 
Bartholome [1], Grilly [36], and Kerr [5] are included with extrap- 
olated values from the data of Bartholome [2] and Friedman et al. 
[6]. The solid curve through these data was calculated from the 
liquid equation of state of Prydz [7]. The lower curve was calcu- 
lated from eq (8). The dashed portion represents an extrapolation 
beyond our measured temperature range. 
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pressure curve. As figure 7 illustrates, these results are 
well described by the Prydz equation. However, 
throughout the plotted range, the earlier values are all 
consistently larger than those extrapolated from this 
work. In the interval 18.7 to 20.4 K, the offset is approx- 
imately 0.5%. About midway between the two curves 
are the extrapolated values from the data of Friedman et 
al. [6] and the direct measurement at 19.48 K by Grilly 
[36]. 

If each research had ±0.15%uncertainty, our results 
agree with [6] and [36] but not with the others. The most 
serious discrepancy appears to be that between Kerr's 
data and ours. Kerr published five measurements of the 
molar volume between 19.5 and 24.2 K. His results at 
19.51 and 21.14 K are larger than corresponding values 
calculated from eq (8) by 0.13 and 0.10 cmVmole, re- 
spectively. Kerr stated that a mass spectrometric anal- 
ysis of the sample showed it to contain 0.4 at % hydro- 
gen (equivalent to 0.8% HD) as the only detectable 
impurity. He corrected the data for impurities in the 
manner described above. If we compare our corrected 
data to his, the differences appear to be as large as 0.12 
to 0.15 cmVmole. No error limits were quoted for the 
results Kerr presented, but if we assume an accuracy of 
±0.15%, these data remain in substantial disagreement 
with our calculated values. 

The source of the large discrepancy is not clear. We 
considered the possibility that some systematic error 
may have arisen from either of the volume calibrations 
at low temperature. Kerr calibrated his pycnometer 
with "pure hydrogen" and the data of Scott and Brick- 
wedde [38], On this point Kerr's report is somewhat 
ambiguous because Scott and Brickwedde included 
measurements on both liquid normal and para hydrogen. 
However, both sets of these earlier data extend along 
the saturated vapor pressure curve from 14 to 20.4 K. 
They were obtained with the same experimental appara- 
tus and, except for sample preparation, by the same 
procedure. We may reasonably assume that they are 
internally consistent. Therefore, regardless of the form 
of hydrogen that Kerr used, we can check his cali- 
bration against our own by comparing the respective 
p-H 2 density standards. 

It was not necessary for us to do this explicitly be- 
cause Roder et al. [39] made this very comparison only 
a short time after the Goodwin data were published. To 
do this, Roder first fitted the Goodwin data to low- 
order polynomial expansions to represent both iso- 
therms and isochores of the liquid density. These 
smoothed functions were then extrapolated to the satur- 
ated liquid-vapor pressure curve reported by Weber et 
al. [40]. Roder finally combined these "derived" data in 
the range 17.0 to 33.0 K with the Scott and Brickwedde 



measurements and fitted the set to a single analytic func- 
tion. They showed that deviations from the fitting func- 
tion were of the same magnitude as the expected error of 
the data. This overall consistency is taken to establish 
the correspondence between our own volume cali- 
bration and that of Kerr's research. Whatever the source 
of the discrepancy is in the deuterium measurements, it 
does not arise from any disparity in the hydrogen densi- 
ties used for the calibrations of the low-temperature 
sample volumes. 

It would appear from figure 7 that Kerr's data agree 
with those of Clusius and Bartholome [1] and Bar- 
tholome [2]. The average deviation of the Clusius and 
Bartholome results from our extrapolated function is 
0.10 cmVmole. But, there are several reasons to suspect 
that these earlier data are too large by roughly 0.05 
cmVmole. The molar volume measurements that Clus- 
ius and Bartholome made were only part of a rather 
comprehensive thermodynamic study. In addition to 
these data and among others, they included mea- 
surements of the triple-point and, at one temperature, 
the saturated liquid-vapor pressure. This information is 
useful because the most likely occurring sample con- 
taminants in the earlier work on deuterium are H 2 and 
HD [33]. It is known that both impurities tend simulta- 
neously to increase measured vapor pressures and to 
depress observed triple-point temperatures. From the 
results of these additional measurements it is argued [33] 
that the sample lot of deuterium used by Clusius and 
Bartholome was contaminated to higher levels than 
were judged in their report. Estimates of either 1.0% H z 
or 2.6% HD are proposed to explain the observed off- 
sets in their results for the vapor pressure and triple- 
point temperature. In either case, the corresponding val- 
ues for AK| calculated from eq (13) are approximately 
0.05 cmVmole. 

In the later high-pressure experiment [2], Bartholome 
calibrated the volume of his cell with both liquid hydro- 
gen and liquid deuterium. Absolute volumes were calcu- 
lated from the measured mole contents by using the 
density data of Onnes and Crommelin [42] for hydrogen 
and Clusius and Bartholome [1] for deuterium. The re- 
sults obtained from the two calibrations reportedly 
agree to within 0.03%. This is interesting because the 
Onnes and Crommelin data are known now to be high 
by about 0.05 cmVmole [38]; the agreement suggests 
that the deuterium results of Clusius and Bartholome 
(and Bartholome [2]) are also high by the same amount. 
The discrepancy is exactly that estimated above from 
considerations of impurity effects. If the Clusius and 
Bartholome data are corrected for this suspected con- 
taminant, their measurements are seen to agree almost 
identically with the direct measurement of Grilly [36] 
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and also with the extrapolated data of Friedman et al. 
[6]. For the Clusius and Bartholome measurements, an 
error limit of ±0.04 cmVmole is not an unreasonable 
estimate, given the combined uncertainties of the data 
and the impurity correction. If this limit is assumed, 
their results are seen to be consistent with our own data 
even after the 0.015 cmVmole impurity correction is 
applied to the latter. 

To summarize this discussion, we have examined 
several series of measurements of the molar volume of 
liquid n-D 2 . Some of these were direct measurements 
along the saturated liquid-vapor pressure curve [1,5,36] 
while others were made along isotherms at pressure 
[2,6]. The individual results were compared directly 
when possible and by extrapolation when not. Gener- 
ally, the values prescribed by our equation of state are 
lower than those obtained from earlier measurements. 
Comparisons with independent measurements of the 
molar volume show our work to be consistent with that 
of Friedman et al. [6], which are the more precise of the 
existing data at pressure. At the vapor-pressure line, we 
find that our extrapolations are consistent with the data 
of Clusius and Bartholome [1] if these are corrected for 
a substantial impurity effect. However, there is a signifi- 
cant difference between the molar volume mea- 
surements of Kerr [5] and corresponding values calcu- 
lated from our equation of state. We examined the 
possibility that some error may have arisen from incon- 
sistencies in the hydrogen density data used as standards 
in the volume calibrations, but we were unable to sub- 
stantiate this hypothesis. The source of the discrepancy 
between Kerr's results and ours is unknown. 

3.6 Other Related Data s 

Most of the previous discussion has involved direct 
comparisons of our data with those of others. There are 
also indirect data comparisons worth considering. One 
quantity that is often discussed is the isochoric tem- 
perature derivative of the pressure, (&P/zT) v . An anal- 
ytic form for (dP/3T) v can be derived in either of , two 
ways: eq (8a) can be inverted to express the pressure as 
a function of temperature and volume, and the desired 
function is obtained by differentiation. Or, the result can 
be calculated directly from the identity 



(dP\ = a 



(14) 



and eqs (9) and (11). In either case, we estimate an 
uncertainty of ±2% for this calculated quantity. 

In figure 8, we have plotted (dP/dT) v for our liquid 
data along the melting curve as a function of tem- 
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Figure 8 — Plot of {t>P/iT) v for liquid n-D 2 along the melting line 
defined by eq (3). The upper solid curve was calculated from eq 
(14), The dashed curve was calculated from eq (15). The lower 
solid curve represents the empirical function for (aP/aT 1 ),, given by 
Prydz [7], 

perature. For comparison, we have also included results 
calculated at the same temperatures and pressures from 
the (dP/dT) p function given by Prydz [7]. The curves 
show two obvious dissimilarities. First, they differ in 
magnitude by approximately 20%. Second, the Prydz 
function shows a large positive curvature, whereas our 
(dP/dT)y results are more nearly constant with only a 
slight positive slope and curvature. The magnitude and 
functional behavior of our (dP/dT)yis consistent with 
that deduced by Friedman etal. [6] from their liquid 
PVT measurements. In particular, they found their data 
best described by the reciprocal volume relationship 



f^ = -14.3 + ^ 



.3 77,- 



Vt 



(15) 



where the molar volume V\ is taken in cmVmole 'and 
(dP/BT)„\s given in bar/K. In figure 8, we have plotted 
the results obtained from eq (15) with the liquid molar 
volume data in table 3. 

Thus far, we have only discussed equation-of-state 
data although these alone are not sufficient to describe 
the thermodynamic properties of a system. Supple- 
mentary thermal data are necessary for complete infor- 
mation. In connection with the present work, it is, there- 
fore, interesting to briefly review the existing 
calorimetric data. 

The earliest specific heat measurements on liquid 
n-D 2 were published by Clusius and Bartholome [1]. 
Their data, over the temperature range 19.4 to 21.7 K, 
were later found to be about 8% lower than those of 
subsequent measurements [35,43]. Kerr et al. [43] con- 
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eluded that some systematic error may have affected the 
earliest liquid results. Today, the liquid heat capacities 
of Clusius and Bartholome remain a matter of historic 
interest but are no longer given as preferred data [44]. 

Subsequent calorimetric measurements were reported 
by Kerr et al. [43] on o-D 2 (2.2% para) in the range 19.9 
to 22.8 K and by Grenier and White [35] on 78.6% para 
at three temperatures between 19.9 and 22.1 K. The 
results of these two experiments agree to within 1 or 
2%. Grenier and White considered this discrepancy rea- 
sonable in view of the uncertainties involved in the inde- 
pendent investigations. They further concluded that any 
effects of ortho-para composition are negligible in this 
temperature and pressure range. 

Still later, Brouwer et al. [45] measured the specific 
heat of liquid n-D 2 as part of an extended study of the 
liquid phase separation properties of the n-Dj/Ne binary 
system. Their results for pure n-D 2 apparently extend 
over the entire range 20 to 30 K because they present a 
difference plot of the raw data from a fitted smooth 
curve. Curiously, Brouwer et al. did not explicitly in- 
clude the analytic function in their publication, although 
they did tabulate values that were calculated from it at 
0.2 K intervals between 24.0 and 27.4 K. 

In each of these four experiments, the actual measured 
quantity is the heat capacity along the saturated vapor 
pressure curve, which we denote as C at . From these 
data, the familiar isobaric quantity C p can be derived 
from the relationship: 



c^rfflffl. 



-'P 



(16) 



where (dP/dT) vp is the slope of the vapor pressure 
curve. The second term in eq (16) may be calculated 
from eqs (8a) and (11) and the temperature derivative of 
Grilly's [23] saturated liquid-vapor pressure function. 
For temperatures below about 23 K, the second term 
in eq (16) contributes less than 1% and could be ne- 
glected, but at higher temperatures there is a significant 
difference between C p and C^, We, therefore, calcu- 
lated C„ explicitly for each of the later sets of mea- 
surements discussed above [35,43,45]. The results were 
then fitted to a quadratic function to yield 

C P (T,P V )=22. 16+0.73(7'- 18.73) 

+0.044(r-18.73) 2 . (17) 

We find that eq (17) reproduces the "corrected" Bro- 
uwer et al. [45] tabulated data to well within 1 %. The C p 
values corresponding to the Kerr et al. [43] mea- 
surements likewise fall in this range except, perhaps, for 



the highest temperature result, which may be lower by 
slightly more than 1%. However, at roughly this same 
temperature, the Grenier and White [35] measurement is 
seen to fall about 2% higher than the smooth curve. We, 
therefore, postulate that eq (17) describes to within 
±1% the isobaric heat capacity for the liquid along the 
saturated vapor pressure curve from the triple point to 
27 K. 

Corresponding values at higher pressures are calcu- 
lated from 

C p (T,P)=C p (T,P vp )- f (d 2 F/dT%dP (18) 

J p v 

where a closed-form expression for the second term is 
derived from eq (8). Equation (18) illustrates that the 
qualitative effect of increasing pressure is to decrease 
Cp . For example, at T = 20.4 K and P = 70 bar, the calcu- 
lated heat capacity is about 9% lower than its value at 
the same temperature on the vapor curve. 

There are no independent measurements of C p at pres- 
sure that could be checked against eq (18) directly, but 
there are data for C y , the specific heat at constant vol- 
ume. Bartholome and Eucken [3] measured C^for liquid 
n-D 2 in the temperature range 19.5 to 23.5 K. In prin- 
ciple, these results could be compared with C sat and our 
equation of state by calculating 



C V (T,V) = C P (TJ 3 )-T V(T,P) 



a\T,P) 
P(T,P) 



(19) 



where the first term on the right-hand-side is obtained 
from eqs (17) and (18) and the second from eqs (8a), (9), 
and (11). 

There is, unfortunately, some ambiguity in the inter- 
pretation of the earlier results. To understand this, it is 
helpful to think of Cp-in terms of temperature and molar 
volume as the independent variables. The difficulty is 
then clear because Bartholome and Eucken did not ex- 
plicitly include the sample density along with their tem- 
perature and specific-heat measurements. From the re- 
ported volume of the sample chamber (3.5 cm 3 ) and the 
amount of material used for the liquid measurements 
(0.171 moles), we might infer the value Fi= 20.47 
cmVmole, but at this density we see from figure 6 that 
a sample could not exist in the liquid state over much of 
the quoted temperature range. We observe from eqs (3) 
and (8) that at 19.65 K the liquid molar volume must 
assume some value between 22.655 cmVmole at the 
melting line and 23.3 10 cmVmole at the saturated liquid- 
vapor pressure. 

Despite this ambiguity, there is some useful informa- 
tion available. The measurements of Bartholome and 
Eucken [3] are plotted in figure 9. If we assume these 
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Bezuglyi and Minyafaev [46] published the first ultra- 
sound velocity measurements on liquid n-D 2 . Their re- 
sults along the vapor curve at T=19 and 20 K are 
plotted in figure 10. Subsequent measurements by Wan- 
ner and Meyer [47] extend from the triple point both 
along the melting line to about 24 K and along the 
saturated-vapor-pressure curve to nearly 26 K. Data 
from the lower portion of their measured temperature 
ranges are included in the plot. Typical error limits for 
these experiments are 0.5 to 1.0%. 

To compare with these direct measurements, we have 
included the results calculated from eq (20) and our own 
data. Sound velocities corresponding to our measured 
temperature range are indicated in figure 10 by the solid 
line segments. The dashed portions represent calculated 



1300- 



Figure 9 — The isochoric heat capacity data of Bartholome and Eucken 
[3], The solid curve was calculated from eq (19) for the liquid 
density at 19.65 K on the melting line. 

data were terminated by the onset of freezing at 19.65 K, 
a sample volume of 22.655 cmVmole is implied. The 
solid line plotted in figure 9 is the C y curve prescribed 
by eq (19) for this molar volume. 

The calculated curve agrees reasonably well with the 
direct measurements. In the range 19.7 to 21.0 K the 
heat capacity results are reproduced to within 5%. At 
higher temperatures the curve follows the mea- 
surements and seems to corroborate the observed weak 
temperature dependence that Bartholome and Eucken 
observed in deuterium. This agreement is somewhat for- 
tuitous, however. The slope and curvature of the func- 
tion defined by eq (19) is very sensitive to the values 
chosen for b 2 and b 3 in eq (1 1), much more so than are the 
results obtained from fitting our data. Below about 20 K 
the calculated C v remains little affected by the choice of 
these two parameters, but at temperatures near 24 K, the 
results can vary by as much as 20%. 

Finally, we compare experimentally measured sound 
velocities in the liquid with those calculated from 



± 2 =m(^1 



V(T,P) 



Ta\T,P) \ 
Cp(J,P) J 



(20) 



where A/=4.0282 gm/mole is the molecular weight of 
deuterium, and u is the longitudinal sound velocity. The 
term in parentheses is calculated from eqs (8a), (9), (11), 
(17), and (18). 
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Figure 10— Liquid sound velocity measurements by Bezuglyi and 
Minyafaev [46], Wanner and Meyer [47], and Pashkov and Kon- 
ovodchenko [48]. Also included in the plot is the semiempirical 
result that [48] derived from multistructure theory and the data of 
Prydz et al. [g]. Equation (20) was evaluated both along the satur- 
ated vapor pressure curve and at melting. The solid lines represent 
these results in the temperature range of the present experiment. 
Extrapolations of eq (20) are shown as the dashed portions. 

values in the extrapolated regions. We estimate roughly 
±2% uncertainties for these results. In the range below 
20.5 K, the experimental measurements and the calcu- 
lated curves agree. There is even a reasonable corre- 
spondence among the results along the vapor curve as 
high as 22 K. However, along the melting line in the 
extrapolated range above 20.5 K there is a significant 
divergence. This discrepancy may indicate a systematic 
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error in our measurements, but it may also be an artifact 
introduced by the particular functional form chosen to 
represent our equation of state. 

Finally, we mention the work of Pashkov and Kon- 
ovodchenko [48] who measured sound velocities at 
three temperatures: 22.00, 26.00, and 29.00 K. At 22.00 
K they obtained 998.8, 1 105.0, and 1 195.8 m/sec at pres- 
sures of 0.605, 50.66, and 101.33 bar, respectively. Cor- 
responding velocities calculated from eq (20) are 983, 
1114, and 1241 m/sec, respectively. At the two lower 
pressures there is satisfactory agreement between the 
measured and calculated values, but at 101.33 bar the 
difference is nearly 4% and is in the same sense as the 
divergence indicated at high temperatures along the 
melting line. 

3.7 Volume Change on Melting Heat of Fusion, 
and Solid Molar Volume 

The volume change on melting (freezing) AV m was 
measured by sealing the sample chamber with liquid 
n-D 2 . At constant temperature, the lower cell chamber 
was first pressurized to freeze the sample and then de- 
pressurized to melt it. The volumes of the solid and 
liquid at the melting pressure were determined quite 
precisely from the abrupt slope discontinuities in P v as 
observed in the calibrated capacitance measurements. 
The results presented in table 4 are expressed as 
APVUm, the volume change normalized to the liquid 
molar volume at the melting pressure. An uncertainty of 
about ±0.0005 is estimated for these data. 

Equation (8) was used to calculate A V m , the absolute 
volume changes. The results are plotted in figure 11 
together with the single value, 2.66 cmVmole, that Clus- 
ius and Bartholome [1] quoted for the triple point. Bar- 
tholome [2] published the only other direct mea- 
surements at these temperatures. His average values, 
2.16 and 1.98 cmVmole at 20.31 and 20.97 K re- 
spectively, are much smaller than our data. The AV m 
measurements by Liebenberg et al. [10] were carried out 
in a much higher range, 75 to 164 K. The results were 

Table 4. Measured volume changes on melting (freezing). 





T m (K) 


A VJ ^(freezing) 
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Figure 11 — Measurements of the volume changes on melting. In- 
cluded with these is the triple-point value derived by Clusius and 
Bartholome [1] from their measurement of the heat of fusion. The 
dashed line is an extrapolation of the empirical function given by 
Liebenberg et al. [10] to represent their high-temperature data. 

fitted to an empirical function of temperature. The 
dashed curve in figure 11 represents this function ex- 
trapolated to the triple-point region. The agreement be- 
tween Liebenberg's extrapolation and the Clusius- 
Bartholome measurement is fortuitious. 

The data presented here are relevant to some inter- 
esting questions first raised more than a decade ago. In 
1973, Roder [49] speculated that hydrogen may exhibit 
a solid phase transformation near the melting curve at 
some temperature between 15 and 26 K. The suggestion 
was made as an attempt to explain certain anomalies 
observed in solid molar volume [50] and heat capacity 
data. Independent measurements of the polarizibility 
[30,51,52] and sound velocity [46,47] in the solid phase 
also showed inconsistencies that seemed to corroborate 
his hypothesis. 

At the same time that Roder's work appeared, Man- 
zhelii et al. [53,54] found that solid hydrogen undergoes 
an abrupt, but reversible, increase in density with in- 
creasing temperature at about 14 K and 30 bar. They 
attributed the effect to a solid phase transition. For a 
time, the proposed explanation drew some criticism 
[13, 14], but supporting evidence was later obtained from 
detailed measurements of the hydrogen melting line 
[24]. A small cusp-like singularity was observed near 14 
K and interpreted as the intersection point of the solid- 
solid phase line. Mills [55] concluded that the equi- 
librium line of the solid hydrogen transformation must 
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have negative slope, but the structure of the high- 
temperature form is still uncertain. X-ray diffraction 
studies at 16 K and 100 bar [56] suggest the change may 
be to an fee structure, but the possibility of an iso- 
morphic hep transformation has not been dismissed en- 
tirely [24]. 

Our volume-change measurements are relevant in this 
context because the solid phase transition is apparently 
not limited to pure p-H 2 alone. Manzhelii et al. [54] and 
Bereznyak and Sheinina [24] discovered that the phe- 
nomenon is independent of the ortho-para concentration 
of hydrogen. Observed irregularities in the melting line 
of e-D 2 [24] further suggest that similar transformations 
may be characteristic of the heavier isotopes as well. 
The effect is apparently also present in hydrogen-helium 
mixtures [57]. It is, therefore, reasonable to expect some 
indication of a solid phase transformation in the melting 
properties of n-D 2 . Yet, as we mentioned earlier, there 
are no anomalies observed in our measurements of the 
melting curve. Nor do we observe any jump-like singu- 
larities in the data for the volume change on melting. 

Manzhelii et al. [53,54] measured a density discon- 
tinuity of 0.15% in solid hydrogen. If a comparable 
difference in density were present between two solid 
phases of n-D 2 , we would expect to see a step in the 
curve representing our results in figure 11. The high- 
temperature solid phase should, by analogy, be more 
dense. The volume change from the liquid would be 
greater, and the step would be positive with increasing 
temperature. The discontinuity would occur at about 19 K 
and be about 0. 3 cmVmole in magnitude.There is no 
obvious indication of such a step discontinuity in figure 
11. Although it could be that the precision and density 
of our data are not sufficient to resolve the effect, it is 
difficult to dismiss the observed regularity of the melt- 
ing pressure measurements. Our result does not directly 
contradict any existing data. At this time there is no 
evidence- for a solid-solid phase transformation in solid 
n-D 2 . 

From measurements of the melting curve and the ac- 
companying volume changes, the heat of fusion AH ( of 
n-D 2 can be calculated with the Clausius-Clapeyron 
equation 






(21) 



Dwyer et al. [58] directly measured the heat of fusion of 
solid p-H 2 and found their results to be well described by 
a linear function of the melting pressure. Although they 
offered no explanation for this simple behavior, it sug- 
gests that we attempt a similar representation. A least- 
squares fit to our calculated heats of fusion gives 



Atf f = 197.22+0.179^ 



(22) 



where the melting pressure P m is in bars and AH f is in 
J/mole. The average absolute deviation of the data from 
eq (22) is 0. 17%. No better fitting results are obtained by 
including quadratic or other high-order terms to the 
empirical function. 

Calculated Aff f are plotted as a function of pressure in 
figure 12. Included in the plot is the earlier direct result 
by Clusius and Bartholome [1], which is seen to be con- 
sistent with the present measurements. An extrapolation 
of eq (22) to the triple point yields a value of 
197.25+0.40 J/mole. Kerr et al. [43] obtained 
196.94+0.40 J/mole for the heat of fusion of o-D 2 at the 
triple point. Measurements on 80% p-D 2 by Grenier and 
White [35] give an average value of 197.41 ±0.40 J/ 
mole. The general agreement among these independent 
data indicates that beyond the nominal error limits of 
+0.40 J/mole there is no measurable dependence of the 
heat of fusion on ortho-para concentration in deuterium. 
The only inconsistent results are those presented by Bar- 
tholome [2] whose data on n-D 2 show Ai/ f decreasing 
with increasing temperature. 

The molar volume V sm of the solid along the melting 
line can be derived from the equation 



Fsm_F|n Y O 



(23) 
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Figure 12— Data for the latent heat of fusion derived from eq (22). 
Included with these is the single direct measurement on n-D 2 by 
Clusius and Bartholome [I]. 



242 



where the liquid molar volume V ]m at the same tem- 
perature and pressure is determined from eq (8). Calcu- 
lated V sm corresponding to the measured values of 
A V m / V im in table 4 are plotted in figure 1 3 . Also included 
in the; plot are two points that Woolley et al. [4] derived 
from the work of Clusius and Bartholome [1] and Bar- 
tholome [2]. 



pUTJ = 0.039 1 33 +0.00053627L. 
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Figure 13 — Data for the molar volume of solid n-Dj along the melting 
curve derived from eq (23). The dashed line was derived from the 
equation of state of Driesseii et al. [14]. Also included are two 
points compiled by Woolley et al. [4] and two neutron diffraction 
results [64,65]. 

The functional behavior of these data suggests a sim- 
ple linear relationship. A least-squares fit gives 



r s m(i > m )=2O.337-O.0053i' ni 



(24) 



where V im (P m ) is in cmVmole and the melting pressure 
P m is in bars. The average absolute deviation of the data 
fronreq (24) is 0.015%. No improvement is obtained by 
including a quadratic term in the fitting polynomial. We 
estimate the accuracy of this^ expression to be about 
,±0.04 cmVmole. : , f 

Within our range of measurements, eq (24) describes 
the data quite well, but it cannot be used reliably for 
extrapolation. A better representation is obtained by fol- 
lowing some of the earlier work on p-H 2 . Dwyer et al. 
[59] first measured the molar volume of solid p-H 2 along 
the melting line at pressures up to 340 bar. Younglove 
[30] later corrected portions of these data and fit the 
solid densities to a linear function of melting tem- 
perature. A similar least-squares fit to the present data 
gives 



where p sm {T) is the density of solid n-D 2 along the melt- 
ing curve in moles/cm 3 . Equations (24) and (25) both 
describe our measurements with comparable accuracy 
and precision, but eq (25) may provide a useful repre- 
sentation for the solid density to temperatures as high as 
27 K. In the following discussion, we examine the con- 
sistency of eq (25), first with independent measurements 
on the low-pressure solid and then with various extrapo- 
lated data at higher pressures. We conclude this section 
with a discussion of the existing measurements of the 
dielectric constant of the solid and liquid phases. 

An extrapolation of eq (25) to the measured triple 
point yields a density corresponding to V im (T tp = 18.7067 
K)= 20.342 ±0.04 cmVmole. It is possible to compare 
this value with the result of existing x-ray [60], electron 
[61], and neutron scattering experiments [62-65] on the 
low-pressure solid. 

For our discussion, the most important of this earlier 
work is the single-crystal neutron diffraction study of 
Nielsen [64]. These data include lattice parameter mea- 
surements on .o-D 2 along its saturated vapor pressure 
curve from 5 to 18 K. From the hep lattice constants 
that Nielsen reported at 5 K, we deduce a molar volume 
of 19.94 cmVmole. With increasing temperature, expan- 
sion of the material in the a and c directions was ob- 
served to be uniform within the limits of experimental 
uncertainty. Nielsen gave an analytical expression for 
the thermal expansion; from it we calculate for o-D 2 a 
molar volume of 20.41 cmVmole at 18.7 K. Next, we 
apply the prescription of Driessen et al. [14] to calculate 
the small volume difference between ortho and normal 
deuterium at this temperature and density. The final 
result for solid n-D 2 , 20.39 cmVmole, is larger by about 
0.05 cmVmole than that found by extrapolating eq (25). 
However, if the lattice measurements are accurate to 
±0.05%, the uncertainty in the corresponding molar 
volume is ±0.03 cmVmole. A small additional error is 
involved with the extrapolation to 18.7 K. When these 
uncertainties are combined with those of the present 
study, the neutron diffraction result is shown to agree 
with our triple-point value. 

Nielsen [64] published the only measurements on 
large single-crystals. The remaining work includes 
powder-diffraction data [60,62,63,65] and one set of 
measurements on relatively thin samples grown on a 
cold substrate [61], We have included a compilation of 
these results in table 5. All of the data were taken at 
temperatures lower than the triple point, but we have 
included corrections to the corresponding solid molar 
volumes of n-D 2 at 18.7 K. in the manner described 
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Table 5. Compilation of existing measurements of the molar volume of solid deuterium. 



References 


Sample 
{% para) 


Temperature 


Measured 

molar volume 

(cmYmole) 


Molar volume 

corrected to 18.7 K a 

(cmVmole) 


Mucker et al. (1965) [62] 
neutron powder diffraction 


33% 
2% 


13 K 
13 K 


20.08 (±0.12)" 
20.08 (±0.12)" 


20.39±0.12 
20.41 ±0.12 


Schuch and Mills (1966) [60] 
x-ray powder diffraction 


33-65%° 


1..5K 


19.80±O.2O d 


20.28+0.20 


Bostanjoglo and Kleinschmidt (1967) 
[61] electron powder diffraction 


33% 


4.2 K 


19.90±0.11 


20. 37 ±0.11 


Mucker et al. (1968) [63] 
neutron powder diffraction 


3% 
33% 
63% 
80% 
80% 


12.9 K 

10.2 K 

2.0 K 

10.5 K 

4.2 K 


20.02 (±0.12)= 
19.91 (±0.12)= 
19.74 (±0.12)= 

19.84 (±0.12) c 

19.85 (±0.12)= 


20.33±0.12 
20.33 ±0.12 
20.32±0.12 
20.37±0.12 
20.46±0.12 


Nielsen (1973) [64] 
single-crystal neutron diffraction 


2% 


5.0 K 


19.94 (±0.03/ 


20.39±0.04 


Yarnell et al. (1975) [65] 
neutron powder diffraction 


2% 


4.2 K 


19.994±0.013 


20.45 ±0.03 



"This calculation is done in three steps. First, the measured molar volume is corrected according to the prescription of Driessen et al. [14] to 
the value corresponding to zero percent para at the same temperature. Next, the temperature correction is made to 1 8.7 K using the analytic 
expression for the linear thermal expansion by Nielsen [64], Finally, the o-D 2 volume is corrected at 1S.7 K to that of n-D 2 again according to 
reference [14], We assume a error of ±0.03 cmVmole for this calculation. The accuracy limit in column 5 is taken as the square root of the sum 
of the squares of the error limit in column 4 and the assumed error in the temperature correction. 

b Estimated error assumes ±0.005 A uncertainty in a and ±0.005 in c/a. 

c No trend was observed in the hep molar volume with para concentration. 

J We set the uncertainty limit equal to the observed data scatter. 

The ±0.6% volume error was calculated from the quoted ±0.2% uncertainty in the nearest neighbor distance. 

The ±0.15% volume error was calculated with an assumed 0.05% uncertainty for the lattice parameter measurements. 



above. Generally, the corrected data are consistent with 
the extrapolated results of the present measurements. 
However, there are relatively large experimental uncer- 
tainties associated with most of the earlier work [60-63]. 
The data of Nielsen [64] and Yarnell et al. [65] are more 
precise. As we mentioned above, the solid molar volume 
obtained by extrapolating eq (25) to the triple point 
agrees with Nielsen's measurements, but it is 0.11 
cmVmole lower than the corresponding value inferred 
from the work of Yarnell. 

At higher pressures, the most instructive comparisons 
are made using the solid equations of state established by 
Driessen et al. [13, 14]. This work was based in large part 
on their measurements of the isochores of solid n-H 2 and 
n-D 2 to pressures up to 2 kbar at the melting line. The 
molar volumes were not measured directly but were 
determined by correlating the isochore data with exist- 
ing measurements at 4.2 K for both H 2 and D 2 [12] and 
along the melting line for H 2 [59]. The equations of state 
for p-H 2 and o-D 2 are presented in tabular form, and 
corrections to arbitrary ortho-para mixtures are accom- 
plished by straightforward calculation. 



A direct comparison of our data to those of Driessen 
et al. is not possible. Their lowest density isochore cor- 
responds to 19.851 cmVmole, which is slightly beyond 
the range of our data. However, they did extrapolate 
their equation of state to low pressures. We interpolated 
their tabulated values to calculate the molar volume of 
n-D 2 at the temperatures and pressures corresponding to 
our melting line. The results are plotted in figure 13. The 
±0.42% uncertainty was taken from figure 22 of their 
paper [14], Our measurements agree with the extrapo- 
lated results of Driessen et al. within the combined error 
limits. However, we note that the most reliable of the 
low-pressure diffraction work [64,65] is more nearly 
consistent with the extrapolated equation-of-state than it 
is with eq (25). 

Driessen et al. quoted a value of 19.95 cmVmole for 
o-D 2 at zero temperature and pressure. This corre- 
sponds closely to the volume, 19.94+ 0.03 cmVmole, 
that Nielsen measured directly at 5 K. According to the 
equation of state, this isochore should intersect the melt- 
ing line at 90.7 bar. The corresponding melting tem- 
perature, 20.90 K, was used in eq (25) to calculate an 
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extrapolated density, and from this we obtain the molar 
volume 19.87±0.04 cmVmole. Nielsen's measurement 
corrected to the normal spin mixture at 20.90 K becomes 
19. 92 ±0.03 cmVmole. Again, the two results agree 
within experimental error. 

In addition to his low-pressure data, Nielsen measured 
the lattice parameters of hep o-D 2 at 5 K and 275 bar. 
From these lattice data, we calculate a molar volume of 
18.64+0.03 cmVmole, which is consistent with the 
Driessen et al. analysis at that temperature. (Existing 
P-V data at low temperatures [12] were fitted to the 
Birch equation. At this solid density, this equation yields 
a pressure of 271.5 bar.) 

According to eq (19) of reference [14], the isochore at 
18.64 cmVmole intersects the melting line at 422.5 bar. 
At this pressure the melting temperature, 27.45 K, in eq 
(25) yields a molar volume of 18.57+0.04 cmVmole, 
which agrees with 18.63+0.03 cmVmole, the Nielsen 
volume corrected to 33% para at 27.5 K. Equation (25) 
is thus seen to reproduce the extrapolated results of the 
neutron diffraction measurements within the combined 
limits of experimental error for the two experiments. 
This equation may, therefore, provide a useful repre- 
sentation for the solid density along the melting line to 
temperatures as high as 27 K. 

There have been a few measurements of the dielectric 
constant of solid [52] and liquid [66,67] n-D 2 . It is in- 
structive to discuss these results in terms of the molecu- 
lar polarizability of the material. To do this, it is neces- 
sary to combine the dielectric data with independent 
measurements of the sample density. The polarizability 
P E can then be calculated from an expression such as the 
Clausius-Mossotti equation 



P _JLJ^v\ 

rE ~ 4ttN\<l+2) 



(26) 



Their value, e= 1.3324, was obtained from a single sam- 
ple under the melting pressure at 21.1 K. The molar 
volume of the sample is calculated from eq (25) to be 
F sm = 19.825 cmVmole. These values substituted into eq 
(26) give P E =7.8388 XlO -25 cm 3 which is about 1.2% 
lower than the values derived from the liquid data. 

The large difference between the derived polar- 
izabilities for the solid and liquid phases is probably 
artificial. The melting pressure at 2 1 . 1 K is only 99.4 bar 
according to eq (3). There should be no significant 
change in the liquid polarizability at this temperature 
over the limited pressure interval between the vapor and 
melting curves. Udovidchenko and Manzhelii [51] mea- 
sured isothermal compressibilities of solid p-H 2 to pres- 
sures as high as 180 bar. They assumed a pressure- 
independent molecular polarizability in their data 
analysis. The results agree with the compressibilities 
that Driessen et al. [14] later found from their isochoric 
measurements. 

Again, there is probably little if any difference be- 
tween the molecular polarizabilities of the solid and liq- 
uid phases at 21. 1 K on the melting line. Once more, our 
justification derives from independent work on p-H 2 . 
Younglove [30] first measured the dielectric constants of 
solid and liquid p-H 2 along the melting line. Although, 
these data suggest a possible difference of 0.4%, sub- 
sequent measurements by Wallace and Meyer [52] and 
Udovidchenko and Manzhelii [51] have not supported 
Younglove's findings. The later experiments obtained 
solid- and liquid-phase polarizabilities that are identical 
to within 0.1 to 0.2%, an accuracy limit imposed by the 
available solid density data. The large apparent differ- 
ence in the deuterium polarizabilities suggests to us that 
more careful measurements of the dielectric constant 
are required for both solid and liquid in this pressure 
range. 



where N i = 6. 0225X10 23 molecules/mole and V is the 
molar volume. 

We first apply eq (26) to the existing liquid data. Ko- 
gan et al. [66] measured the dielectric constant of liquid 
n-D 2 along the saturated vapor-pressure curve from the 
triple point to 20.4 K. From their 12 data points and the 
corresponding liquid molar volumes computed from eq 
(8), we calculate an average value of P s =7.9371 X 10 -25 
cm 3 with a standard deviation of 0.0014X 10 -23 cm 3 . 
This result agrees with those derived from similar data 
reported by Constable et al. [67] in the range 18.9 to 23.0 
K. From their 22 measurements, we calculate an aver- 
age of 7.9330 X 10 _2S cm 3 with a standard deviation of 
0.0077xl0- 25 cm 3 . 

To our knowledge, the only solid phase measurement 
on n-D 2 is that published by Wallace and Meyer [52]. 



3.8 Solid Thermal Expansivity 
and Isothermal Compressibility 

Measurement accuracies of solid-phase properties are 
severely constrained with the present technique. The 
problem derives from the nonplastic-flow character- 
istics of the solid and a rather sensitive dependence of 
the experimental results on the homogeneity of the sam- 
ple density. For measurements of both thermal expan- 
sivity and isothermal compressibility, it is necessary to 
measure changes in the sample volume that result from 
unequal deflections of the upper and middle diaphragms 
of the pressure cell (see Fig. 1). Volume changes of the 
thin disk-shaped solid always introduce radial variations 
in the compression. As P\ is increased, for example, the 
middle diaphragm is displaced upward, and the sample 
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is compressed to a higher density at the center than 
around the periphery. 

For liquids, the distortion of the sample geometry is 
no problem. Reasonably accurate results were even ob- 
tained earlier on solid helium [16,17], although hyster- 
esis effects, particularly at the lowest temperatures indi- 
cate similar problems. With solid n-D 2 , we cannot 
assume that our samples were free of internal strain even 
after several hours equilibration near the two-phase re- 
gion. We observed that density gradients persisted in 
sufficient degree to preclude all but very approximate 
measurements of a and /3 for the solid. 

We tried a number of approaches to defeat the prob- 
lem but could only obtain reproducible results by allow- 
ing long (2-10 h) relaxation periods with the solid at 
temperatures and pressures near the melting curve. Typ- 
ically, in the melting pressure range 20 to 70 bar, the 
isothermal compressibilities decreased from about 3.8 to 
2.0X10 -4 bar -1 with increasing P m . Thermal expan- 
sivities were in the range 2.4 to 2.6X 10~ 3 K _1 for melt- 
ing temperatures between 19.2 and 20.4 K. 

These results are somewhat inconsistent with our 
solid molar volume data. Consider the equation 

V sm \dP m J dP m /dT m 

where & m and ct sm denote respectively the isothermal 
compressibility and thermal expansivity for the solid at 
melting. The quantity on the left-hand side of this equa- 
tion represents a coefficient of compression for the solid 
along the melting curve. It can be calculated from the 
fitted curve of eq (24); the resulting values increase from 
2.61 X 10 4 bar - ' at the triple point to 2.65 X 10~ 4 bar -1 
at 70 bar. The average measured values over the pres- 
sure range 20 to 70 bar, /3 sra =2.73xl0^ 4 bar -1 , 
a sm = 2.56x 10- 3 K-\ and d/ , D /dr ni =41.5 bar K" 1 , yield 
the value 2.1 1 X 10 -4 bar -1 for the right hand side of eq 
(27), 

Estimates of the solid-phase properties are available 
from the Driessen et al. [14] equation of state. Although 
their tabulated equation is for o-D 2 , the compressibility 
and thermal expansion coefficient for n-D 2 can be calcu- 
lated from the information provided. In the temperature 
range 18.7 to 21.0 K, the differences between the differ- 
ential quantities of o-D 2 and n-D 2 are less than 1%. At 
the triple point, we calculate a compressibility for solid 
n : D 2 of 4.28 Xl0~ 4 bar^ 1 and a thermal expansion of 
about 5.3X 10~ 3 K" 1 . At T m =20.5 K, the corresponding 
values are about 12% lower. 

An indirect estimate of the solid thermal expansivity 
can be made from some of the data that we have already 
discussed. Consider the temperature change of the liq- 
uid enthalpy along the melting curve 



dffim r v dP, 
&T~ PM ^ lm dr n 



'■n ( \ 

r 1 [ 1 — ?V*lml- 



(28) 



In this expression C PM denotes the isobaric heat capacity 
of the liquid at the melting line. The right-hand side of 
eq (28) can be evaluated from eqs (4), (8), (11), and (17). 
At the triple point we obtain d^i ra /d7* m = 9 1.5 ±1.0 J/ 
mole K. 

The enthalpy of the solid at melting H sm is equal to the 
liquid enthalpy minus the heat of fusion. We may, there- 
fore, write 



dH im _dH lm _ d(^H ( )dP m 
dT m dT m di> dr m 



(29) 



The second term evaluated from eqs (4) and (22) con- 
tributes an amount 7.0+0.5 J/mol K. The combined 
results in eq (29) gives dff an /dr m =84.5±1.2 J/mol K. 

Clusius and Bartholome [1] and Kerr et al. [43] mea- 
sured the specific heat of the solid along the vapor pres- ' 
sure curve. These data extrapolate to the triple point to 
yield C PiSm — 11.0+0.1 J/mole K. This value combined 
with dH sm /dT m , V im from eq (25), and dP m /dT m at the 
triple point give a sra =3.8±0.6x 10~ 3 K -1 . This indirect 
estimate of the solid thermal expansivity lies about mid- 
way between the value derived from the work of 
Driessen et al. [14] and our direct result. 

We began this section with a brief discussion of the 
problems that we met in applying our technique to solid 
samples. It is appropriate to conclude with some sug- 
gestions for improving the method. As we mentioned 
earlier, the major problem derives from the design of 
our pressure cell. Impressed changes in the sample vol- 
ume always introduce radial inhomogeneities in the den- 
sity. For solid deuterium, the very thin disk-shaped ge- 
ometry of the sample combines with an inherently large 
internal friction. Together, they allow substantial pres- 
sure gradients to persist for prohibitively long periods of 
time. 

We were unable to improve our results by annealing, 
even when this was done at pressures and temperatures 
very close to melting. We feel the only possible solution 
is to redesign the cell in a way that minimizes the radial 
density gradients. For example, it may be possible to 
construct a cell with a sample chamber thickness that is 
very nearly proportional to the displacement of the di- 
aphragms. If this can be done, it is clear that changes in 
the sample volume will induce density gradients in the 
axial rather than the radial direction. Sample distortions 
can never be totally eliminated, but they should be far 
less severe in a carefully designed, lens-shaped sample 
than they were shown to be with our geometry. 
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3.9 Further Observations 

Finally, we discuss the effects of para to ortho con- 
version and their possible influence on the PVT results. 
We ran a number of experiments to measure such effects 
explicitly in both hydrogen and deuterium. For both 
isotopes, spin compositions were determined from the 
thermal conductivity of the gas samples [68]. On solid- 
liquid mixtures of n-H 2 , the measured time rate of 
change in x p> the para mole fraction, was djCp/dr=0.01 
h _I . This result agrees with the self-conversion rate 
given by Woolley et al. [4]. Measurements of the melting 
pressure P m and the liquid molar volume V\ (close to the 
melting curve) as functions of time gave dP m /dr=0.06 
bar h -1 and dfyp;df=7x 10 -5 h _1 . These results di- 



vided by 



dxp/dr 



gives d/' ra /dXp=6 bar and 



dV\/V] cbc p =7 X 10 i , which are consistent with previous 
results: / , m (p-H 2 )~P m (n-H z ) from the data of Woolley et 
al. [4], Mills and Grilly [28], and Goodwin and Roder 
[31]; Fi(p-H 2 )-F,(n-H 2 ) by Scott and Brickwedde [38] 
and Wallace and Meyer [69]. 

For liquid n-D 2 , the rate of self conversion, djc D /dr, 
was estimated by Woolley et al. [4] to be 5X 10~ s h _1 . 
For the solid, Motizuki [70] calculated the value 5 X 10 -4 
h _1 which is in the range 2 to 8X 10" 4 h _1 that Grenier 
and White [35] obtained by direct measurement. We 
measured the conversion rate of n-D 2 several times in 
our cell. The results for dx /dr were between 1 and 
3X 10~ 4 h _1 for observation periods of 60 to 90 h. The 
data are reasonably consistent with the earlier the- 
oretical and experimental work. The conversion rates 
are small, and we feel that our PVT measurements are 
not significantly influenced. The following numerical 
estimates should illustrate this point. 

One effect of self-conversion is an increase in P m at 
constant temperature. We noted in previous discussion 
that the melting pressures reported by Bereznyak and 
Sheinina [24] on o-D 2 exceed ours on n-D 2 (33.3% ortho) 
by about 2 bar. This difference scales to 6 bar for 
A* p = 1, which is similar to the difference in P m between 
pure ortho and para hydrogen. For the self-conversion 
rate dx„/d*=2X 10~ 4 h _1 , we should observe 
di> m /dr«1.2x 10- 3 bar/h. 

Another effect of conversion is an increase in the 
molar volume of the sample. An estimate of the frac- 
tional change in the liquid volume can be calculated 
from the equation 



v! A U£> 



Ax„ 



(30) 



fraction of thepara spin modification in the mixture, and 
P Q is an effective pressure introduced by the electric 
quadrupole interactions of thepara species. We estimate 
the value (dPq/dx ? ) r ,y^l2.3 bar from table 1 of Driessen 
et al. [14]. For liquid compressibilities on the order of 
7X UT 4 bar -1 , we calculate AJVFi=8x 10 3 for Ax p =\, 
which is also similar to H 2 . Thus in n-D 2 with 
dXp/dt = — dxj&p, we expect the rate 
dH/Fidr«2xlO- 6 h-'. 

The longest times taken for the various measurements 
were 1 h for the liquid a and jS, 10 h for solid a. and /?, 
and 20 h for A V m . The longest time of low-temperature 
exposure before P m and V\ measurements were made was 
12 h. Thus, in all cases the effects of ortho-para con- 
version are within the quoted accuracy limits of the 
data. It is nonetheless interesting to compare the values 
of P m and Vi at the beginning and at the end of the long 
runs (60 to 90 h) involving AF m . 

Again, with an assumed conversion rate of 2X 10 -4 
h" 1 , we expect a fractional volume change in the liquid 
to be roughly 1.6 XlO -4 for an 80-h period. Volume 
changes of this magnitude are easily measured in our 
apparatus. Indeed the observations include such results, 
but there were also some changes that are much smaller 
and, in two cases, some large volume decreases (1.7 to 
1.9X 10 -4 ). However, volume losses were also observed 
in the p-H 2 measurements, where no conversion occurs, 
and these can be explained by a loss of the sample 
through a leaky cell valve. 

For the melting pressure this self-conversion rate 
should produce an increase of about 0.10 bar for n-D 2 
during an 80-h period. Actual measured values, shown 
in table 6, vary from —0.39 to +0.38 bar. The large 
decreases in P m are not easily explained. If they were to 



Table 6. Observed changes in the melting pressures of samples that had 
been held at low temperatures for 50- to 70-h periods. 



r(K) 



Seq. No. 1 



A/ybar) 



18.83 
19.00 
19.20 
19.20 
19.40 
19.60 
19.80 
19.80 
20.00 
20.20 
20.20 
20.40 
20.40 
20.40 



9 

5 

8 

11 

4 
6 
3 

12 

7 

1 

2 

10 

13 

14 



0.00 

-0.39 

-0.28 

0.10 

0.31 

-0.01 

-0.06 

0.19 

-0.30 

0.24 

-0.03 

-0.33 

0.04 

0.38 



Where ft is the Compressibility Of the liquid, X p is the 'Chronological order of the measurements. 
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be ascribed to bath temperature decreases, the A7 1 
would have to be on the order of 10 mK, which is far 
greater than our expected precision of a millikelvin or 
two. 

There is evidence in some cases that a high initial 
measurement of P m could have been made before the 
sample was adequately cooled; however, in other cases 
this explanation is unreasonable. For example, after the 
first run at 19.2 Kwas completed, we saw that a large 
decrease (AP m = —0.28 bar) had occurred in the melting 
pressure. We replaced the original sample with a fresh 
charge of n-D 2 from the storage bottle. The melting 
pressure of the second sample was then measured, and 
the result agreed with that obtained from the first sam- 
ple when it was originally admitted to the cell and pres- 
surized to cause freezing. The result removes the possi- 
bility of some kind of slow cooling of the cell, but it 
introduces the possibility of inadequate time for the liq- 
uid to be cooled or changed in some way with time. 

Despite the lack of evidence for two solid structures 
in n-D 2 , we did consider this possibility in connection 
with our observed melting pressure changes. If the 
Gibbs-energy functions of two solid phases are approxi- 
mately equal, it should be quite easy to pressurize a 
liquid sample to the metastable solid form. It can be 
shown in this case that the observed P m will be higher 
than that corresponding to the stable phase. At constant 
temperature the solid may then revert to the equilibrium 
structure, and as it does, the observed melting pressure 
will decrease. This argument is seen to hold for tem- 
peratures both above and below the assumed solid-solid- 
liquid triple point. 

Although, as we mentioned, there is no direct evi- 
dence for more than a single solid structure in n-D 2 , we 
did observe some curious behavior with p-H 2 in our cell. 
Recall that Manzhelii [53,54] reported volume decreases 
of 0.15% in solid p-H 2 with small increases of tem- 
perature. In p-H 2 and n-H 2 , Bereznyak and Sheinina [24] 
observed kinks in the freezing curve near 14.9 K. Simi- 
larly, in our PVT apparatus, we observed A V=±— 0.07% 
in solid p-H 2 with slight increases of temperature. We 
also found the melting pressure of p-H 2 to decrease with 
time at rates of about 0.005 bar h _1 . Over time periods of 
60 h, P m decreases of approximately 0.3 bar were ob- 
served, similar to our measurements on n-D 2 . 

If this similarity supports the possibility that two (or 
more) solid phases exist in n-D 2 , there is contrary evi- 
dence as well. In one of the long runs at 20.4 K, a 
melting pressure change of —0.33 bar was observed. 
After the run was completed but before the low- 
temperature valve was opened, the sample was held in 
the liquid phase for over an hour to assure that all of the 
solid had melted. The liquid was then refrozen and the 



melting pressure was again measured. The resulting 
value was identical to that obtained from the melting 
process only an hour before. This experiment tends to 
eliminate the possibility of a solid structure change with 
time. At best, it could indicate that such a change is fast 
when an "old" liquid is used. At this time, the anomalous 
changes of the melting pressure are not well understood. 



4. Summary and Conclusions 

The purpose of this study has been to provide an 
accurate and precise determination of the equation of 
state of n-D 2 in the region of the phase diagram near the 
triple point. We began by presenting measurements of 
the melting line. These data range in temperature from 
the triple point to 20.4 K. The melting pressures were fit 
to a quadratic function of temperature. The average 
deviation of the points from the smooth curve is about 
0.1 bar, which is comparable to the precision of the 
measurements. The observed regularity of these data 
contrasts to the anomalous behavior reported in similar 
studies [24] on p-H 2 , n-H 2 , and o-D 2 . The results of the 
earlier work were interpreted as evidence for the exis- 
tence of solid-solid-liquid triple points in these systems. 
While this may be true, there appears to be no firm 
indication of a similar phenomenon in n-D 2 . 

We combined our melting data with existing vapor- 
pressure measurements to establish an accurate triple- 
point temperature. The result was corrected for the HD 
impurity measured by mass-spectrometric analysis. The 
value T Xv = 18.723±0.010 K was given as the triple-point 
temperature of pure-phase n-D 2 . A detailed comparison 
of this result to existing data is given elsewhere [33]. 

Data for the isothermal compressibility and thermal 
expansivity of the liquid were collected in the tem- 
perature range 18.8 to 21.0 K and at pressures between 
4 and 70 bar. We compared these with estimates of 
corresponding values derived from existing molar vol- 
ume data [2,5,6], The most accurate and precise of these 
earlier data [5,6] gave results that agree with ours to 
within a few percent. We also presented measurements 
of the molar volume of the liquid, which were made at 
nine points along the melting line. The combined liquid 
data were then fitted to a single empirical equation of 
state. The purpose of doing this is twofold. First, a de- 
scription of the data by a single function is sufficient 
demonstration of their internal consistency. Second, the 
equation of state can be used to extrapolate beyond the 
limited range of the data and also to predict certain 
quantities that were not measured directly. This allows 
us to make indirect data comparison that would not 
otherwise be possible. 
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As an example of this, we began by interpolating our 
equation of state to compare the molar volume with 
existing measurements [2,6] at pressure. We then extrap- 
olated the equation to the vapor-pressure curve where a 
number of direct measurements of the liquid molar vol- 
ume [1,5] exist. The analysis was finally extended by 
incorporating some of the existing calorimetric data for 
the liquid [35,43,45]. This additional information al- 
lowed us to compare our results with existing sound 
velocity [46-48] and isochoric heat capacity mea- 
surements [3]. The overall conclusion of the analysis is 
that our measurements agree with the best available data 
for differential quantities, such as the compressibility, 
\(dP/dT) v , and sound velocity. However, small but sig- 
nificant uncertainties remain in the scale of the liquid 
data. Our volume measurements agree with those of 
Friedman et al. [6], but extrapolations of eq (8) to the 
vapor curve yield values that are at least 0. 10 cmVmole 
lower than the measurements of Kerr [5], 

Results for the volume change on melting were then 
given. These were followed by a brief review of the 
existing evidence for a solid-solid phase transition in 
hydrogen. The behavior of our volume-change mea- 
surements showed no irregularities that would indicate 
a similar behavior in n-D 2 . 

From the volume-change measurements and the slope 
of the melting curve we derived the corresponding 
heats of fusion and fitted them to a linear function of 
pressure. These results agree with existing mea- 
surements at low pressures [35,43], Finally, we derived 
solid molar volumes along the melting curve. The corre- 
sponding solid densities were fit to a linear function of 
the melting temperature. This empirical function was 
extrapolated to both low and high temperatures to com- 
pare with existing data. We found that our extrapo- 
lations are consistent with values derived from the best 
existing solid-phase data [14,64]. 

We attempted to measure compressibilities and ther- 
mal expansion coefficients in the solid. The results of 
these measurements were found to be somewhat incon- 
sistent with our solid densities along the melting curve. 
We briefly discussed the compressibility and thermal 
expansion coefficients that can be derived from the tabu- 
lated equation of state that Driessen et al. [14] presented. 
These results appear to be more nearly consistent with 
the observed solid compression. Finally, we derived an 
indirect estimate of the solid thermal expansion from a 
variety of data that had been presented earlier in the 
discussion. This calculation gave a higher value than we 
obtained from our direct measurements but lower than 
that deduced from the extrapolated equation of state. 
The paper concluded with our unexplained obser- 
vations of the changes of the melting pressure of the 



confined sample with time. We noted that these changes 
were opposite from what might be expected of the ef- 
fects of ortho-para conversion. 
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